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Abstract 

We review the non-perturbative theoretical framework set up recently to compute 
the inelastic scattering cross section from quantum impurities [G. Zarand et al, 
Phys. Rev. Lett. 93, 107204 (2004)] and show how it can be applied to a number of 
quantum impurity models. We first use this method for the S = 1/2 single-channel 
Kondo model and the Anderson model. In both large plateau is found 

in the inelastic scattering rate for incoming energies above TV, and a quasi- linear 
regime appears in the energy range 0.05 T% < uo < 0.5 Tk, in agreement with 
the experimental observations. We also present results for the 2-channel Kondo 
model, the prototype of all non-Fermi liquid models, and show that there half of 
the scattering remains inelastic even at the Fermi energy. 
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1 Introduction 



One of the major ingredients of mesoscopic physics is quantum interference: it 
leads to phenomena such as weak localization, Aharonov-Bohm interference, 
universal conductance fluctuations, or mesoscopic local density of states fluc- 
tuations [1]. All these phenomena rely on the phase coherence of the conduc- 
tion electrons. This phase coherence is, however, destroyed through inelastic 
scattering processes, where an excitation is created in the environment. These 
inelastic processes suppress quantum interference after the so-called dephasing 
time, Tip, also called inelastic scattering time. The excitations created in course 
of an inelastic scattering process may be phonons, magnons, electromagnetic 
radiation, or simply electron-hole excitations. 

A few years ago Mohanty and Webb measured the dephasing time t v (T) care- 
fully down to very low temperatures through weak localization experiments, 
and reported a surprising saturation of it at the lowest temperatures [2] . These 
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experiments gave rise to many theoretical speculations: intrinsic dephasing 
due to electron-electron interaction [3] as well as scattering from two-level 
systems [4,5] have been proposed to explain the observed saturation, and in- 
duced rather violent discussions [3,6,7]. Recently, it has been finally proposed 
that an apparent saturation could well be explained by inelastic scattering 
from magnetic impurities [8,9]. 

Triggered by these results of Mohanty and Webb, a number of experimen- 
tal groups also revisited the problem of inelastic scattering and dephasing 
in quantum wires and disordered metals: A series of experiments have been 
performed to study the non-equilibrium relaxation of the energy distribution 
function in short quantum wires [10]. These energy relaxation experiments 
could be well explained in terms of the orthodox theory of electron-electron 
interaction in one-dimensional wires [11], and/or inelastic scattering mediated 
by magnetic impurities [12,13,14,15]. Parallel to, and partially triggered by 
these experiments, a systematic study of the inelastic scattering from mag- 
netic impurities has also been carried out recently, where inelastic scattering 
from magnetic impurities at energies down to well below the Kondo scale has 
also been studied [16,17,18]. 

Theoretically, this strong coupling regime can be reached only through a non- 
perturbative approach. Such a method to compute the inelastic scattering 
cross-section has been proposed in Ref. [9] and further developed in Ref [19], 
where it has been shown that the finite temperature version of the formula 
introduced in Ref. [9] describes indeed the dephasing rate that appears in 
the expression of weak localization in the limit of small concentrations. Ex- 
cept for very low temperatures, where a small residual inelastic scattering 
is observed [17,18], these calculations were in excellent agreement with the 
experiments, and they clearly showed that magnetic impurities in concentra- 
tion as small as lppm can already induce substantial inelastic scattering. We 
have to emphasize though that experiments on very dirty metals, e.g., prob- 
ably cannot be explained in terms of magnetic scattering, and possibly other 
mechanisms are needed to account for the dephasing observed at very low 
temperatures in these systems [20]. 

Here we review the theory of Ref. [9] and show how it can be applied to various 
quantum impurity problems. In Ref. [9] we formulated the problem of inelastic 
scattering in terms of the many-body ^-matrix defined through the overlap of 
incoming and outgoing scattering states: 



out(/IOin= m(/l^)in- i 1 ) 

The incoming and outgoing scattering states, |i)j n and |/) ut' are asymp- 
totically free, however, they may contain many excitations, i.e. they are true 
many-body states. The many-body T-matrix is defined as the 'scattering part' 
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of the S-matrix, S — I + if, with I the identity operator. By energy conser- 
vation 



m 



(f\m in = 2nS(E f -E t )(f\T\ i ), (2) 



where we introduced the on-shell T- matrix (f\T\i). The results of Ref. [9] rely 
on the simple observation, that (pcr|T|pV) determine both the total (er tot ) 
and the elastic (<r el ) scattering cross sections of the conduction electrons (or 
holes) at T = temperature. The total scattering cross section of an electron 
of momentum p and spin a is given by the optical theorem as 

"tow = — Im(p<7|T|p<r> , (3) 



where vf denotes the Fermi velocity. This expression accounts also for pro- 
cesses where a single electron scatters into many excited electron states (see 
Fig. 1). In case of elastic scattering, on the other hand, an incoming single elec- 
tron state is scattered into an outgoing single electron state, without inducing 
any spin or electron- hole excitation of the environment. The corresponding 
cross section can be expressed as 

•s^/^^-fli^' (4) 



with £ the energy of the electron measured from the Fermi surface. Inelas- 
tic scattering processes, schematically illustrated in Fig. 1, can be defined as 
scattering processes, which are not elastic. Accordingly, the inelastic scatter- 
ing cross section associated with these processes is just the difference of these 
two cross-sections: 

"inel = "total ~ °~cl ■ ( 5 ) 



This simple formula allows us to compute the inelastic scattering cross-section 
in detail. 

In a free electron gas it is convenient to introduce angular momentum channels, 
L = (l,m), and define the scattering states in terms of radially propagating 
states |p, a) — > ||p|, L, a). In this basis the on-shell S and T-matrices become 
matrices in the quantum numbers L, and they depend only on the energy uo 
of the incoming particle [21], 

SL,a; Lv(w) = $L,L'$cr,cr' + « *L,<r; LV M • (6) 
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Fig. 1. Sketch of (a) elastic and (b) inelastic scattering processes. In case of an in- 
elastic scattering the outgoing electron leaves spin- and/or electron-hole excitations 
behind. 

By unitarity, the eigenvalues s\ of the matrix sl,u-l'u' must all be within 
the complex unit circle for any u, and they are directly related to the in- 
elastic scattering cross section. In case of s-scattering and spin-conservation, 
e.g., SL,a;L'a' becomes a simple number, s(u) — 1+i t(u), and the inelastic 
scattering cross section can be expressed as 

OfcdM = 4p (1 - |sMD = 4" ( 2 Im *M ~ l*( W )| 2 ) > ( 7 ) 

Pf Pf 

where we assumed free electrons of dispersion £ = p 2 /2m — // with a Fermi 
energy /i and a corresponding Fermi momentum p^. Eq. (7) implies that the 
scattering becomes totally elastic whenever s(u) is on the unit circle, and it is 
maximally inelastic if the corresponding single particle matrix element of the 
S'-matrix vanishes. The former situation occurs at Fermi liquid fixed points, 
while the latter case is realized, e.g., in case of the two-channel or the two- 
impurity Kondo models. The total scattering cross section, on the other hand, 
is related to the real part of s(uj) as 

a tot {co) = — (1 - Re{s(u,)}) = 1m{t(u)} . (8) 
Pf Pf 

It is easy to generalize this result to the case of many scattering channels, and 
one finds that inelastic scattering can take place only if some of the eigenvalues 
of Sl i<t . l'o' ar e not on the unit circle [31]. 

To determine t LaL i a i we need to compute the matrix element (per|T|pV), 
that we first relate to the conduction electrons' Green function through the 
so-called reduction formula [21,22], 

(p,a|T|pV> = - [G°]^ >±CT (0 G ±P±CT)±P , ±(J ,(0 [G°] ± l /± AO ■ (9) 

Here the ± signs correspond to electron and hole states with £ > and £ < 
and of energy E = |£|, G denotes the free electron Green's function, and G the 



full many-body time-ordered electron Green's function. By Eq. (9) the positive 
frequency part of the Green's function describes the scattering of electrons, 
while the negative frequency part that of holes. In case of a degenerate vacuum 
state one must average over the various vacuum states in Eq. (9) [31,19]. 

According to Eqs. (3), (4), (5), and (9), to compute the inelastic and elas- 
tic scattering cross-sections, we only need to evaluate the self-energy of the 
conduction electron's Green function, many cases referred to as the T-matrix. 
This can be done either analytically using, e.g., perturbative methods, or nu- 
merically, by relating the self-energy to some local correlation function, and 
computing the latter by Wilson's numerical renormalization group (NRG) 
method [34]. The latter approach enables us to compute both the imaginary 
and real parts of the matrix t(u), and we can thus also determine the complex 
eigenvalue s(cu). 



2 Inelastic scattering in the Kondo model 

Let us first apply this formalism to study the inelastic scattering in the single- 
channel Kondo model defined by the Hamiltonian 

h = Y,£p a L a P« + i $ 2 Op^'V,*' • (io) 
p> ct p,p' 

Here a pa creates a conduction electron with momentum p, spin a, and S — 1/2 
is the impurity spin. The T-matrix of the Kondo model can be related to the 

— * 

Green's function of the so-called composite Fermion operator, F a = J2a',p S ■ 
a aa 'a pa .' [30], whose spectral function can then be computed using NRG [34]. 

Before presenting our numerical results, let us briefly discuss what we can 
learn about the inelastic scattering cross-section from analytical approaches. 
The uj ^> Tk regime is accessible by perturbation theory. Summing up the 
leading logarithmic diagrams we find 

t(u » T K ) a: % (n 2 /2)S(S + l)/\n 2 (iu/T K ) , (11) 

where Tk ~ Epe" l ^ Je is the Kondo temperature, with Ep the Fermi energy 
and q the density of states at the Fermi energy for one spin direction [29] . To 
leading logarithmic order, the scattering is completely inelastic 

a incl (cu > T K ) ^ a tot (iu) ^ \ S(S + 1) 2 1 , (12) 

Pf In (u/T K ) 
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Fig. 2. Left: Renormalization group flow of the eigenvalue of the single particle 
5-matrix for the single-channel Kondo model. Right: Energy-dependence of the 
elastic-, total- and inelastic scattering rates for the Kondo model. 



since the elastic contribution only increases as cr cl ~ |t(u;)| 2 , 

1 



aJu > TV 



TP 



Ap 2 F 



s 2 (s + iy 



\n 4 (u/T K ) 



(13) 



This very surprising result contradicts conventional wisdom, which tries to 
associate inelastic scattering with spin-flip scattering. It can be explained in 
the following way [24]: At high energies, incoming electrons are scattered by 
the impurity spin fluctuations. These fluctuations can absorb an energy of 
the order of ~ TV, and therefore the energy of the incoming electron is not 
conserved even in leading order, but it typically changes by a tiny amount, 
5uj ~ Tk- In earlier approaches, this tiny energy transfer has been neglected, 
and the spin-diagonal scattering has been incorrectly identified as an elastic 
process. 

We can also relate the cross sections above to scattering rates. Assuming 
a finite but small concentration n; mp of magnetic impurities, the conduction 
electrons' lifetime can be expressed as \/t{uj) = n imp v F cr imp (u;). This relation 
can be used to define the inelastic scattering rate too as 

1 / \ 7lS(S + l) , 

= n imp v F (Ti ne i{u) w n imp - — 2 , {u > T K ). (14) 



Tinei 2 g In (lo /Tk) 



Note that the latter asymptotic expression is a factor 3/2 larger than the 
Nagaoka-Suhl formula, which takes into account only spin- flip processes [33]. 

For energies |u;| < T# perturbation theory breaks down, and it is more appro- 
priate to use Nozieres' Fermi liquid theory, according to which scattering at 
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Fermi energy scattering is completely elastic, and is described through sim- 
ple phase shifts [23], t a (uj = + ) = 2smS a e %5a . Here we allowed for different 
phase shifts in the spin up and spin down channels. This Fermi liquid expres- 
sion yields 

0tat,<r(w 0) = — sin 2 (5 CT ) , 0wbI,*( u 0) = . (15) 

Pf 

The maximum total scattering cross section is reached in the unitary limit, 
5 a = n / 2, while the inelastic scattering cross-section always vanishes at the 
Fermi energy. Perturbation theory around the Fermi liquid fixed point predicts 
0) oc {oo/T K ) 2 [9,25]. 

The analytical calculations above can only capture the physics in the limit of 
very large and very small frequencies, and for energies uj ~ Tk we need to use 
more sophisticated methods such as NRG. In Fig. 2 we show the evolution of 
the eigenvalue of the s(u). In the limit T^jEp <C 1 this becomes a universal 
function, s(iv) = s(oj/Tk)- In the single-channel Kondo model — > 1 
for both very large and very small frequencies, and thus scattering becomes 
completely elastic both limits. The reasons are different: At large energies 
conduction electrons do not interact with the impurity spin efficiently At 
very small energies, on the other hand, the impurity's spin is screened and 
disappears from the problem [23]. The maximum inelastic scattering is reached 
when the eigenvalue s(iv) is closest to the origin, i.e., at energies in the range 
of the Kondo temperature, uj ~ Tk- 

The total, elastic and inelastic scattering cross sections of an electron are 
shown in Fig. 3. As expected, the inelastic amplitude always vanishes at the 
Fermi level, and at energies well above T K both the elastic and the inelastic 
scattering cross-sections follow the the analytical expressions, Eqs. (12) and 
(13). For energies <C T K we recover the quadratically vanishing inelastic 
rate expected from Fermi liquid theory [25], but the cr inc i ~ u 2 regime appears 
only at energies well below the Kondo temperature, uj < 0.05 Tk (see Fig. 3). 
For | a; | T K the usual Nagaoka-Suhl expression describes the inelastic scat- 
tering rather well apart from the incorrect overall pre-factor 3/2 discussed 
before, but it starts to deviate strongly from the numerically exact curve at 
approximately 10 Tk, and it completely fails below the Kondo temperature 
T K . 

Interesting features appear also at intermediate temperatures. The inelastic 
scattering rate is roughly linear between 0.05 T K < uo < 0.5 T K) and a broad 
plateau appears above the Kondo scale, where the energy-dependence of the 
inelastic scattering rate turns out to be extremely weak. Even though our cal- 
culation is done at T = temperature, a- me \(T, uj = 0) is expected to behave 
very similarly to <7 ine i(T = 0,uj). Thus both features are perfectly consistent 
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Fig. 3. Inelastic, elastic, and total scattering rates for the single-channel Kondo 
model in units of oo = 4-7T /pp, as a function of the incoming electron's energy. 

with several experiments [26,16,8], and, somewhat surprisingly, our results fit 
the experimentally measured temperature-dependence of l/r$ excellently [16]. 
This is beyond expectations, since realistic magnetic impurities have a com- 
plicated d- level structure, and our T = temperature results provide just 
approximations for the dephasing rate, and one should use the finite temper- 
ature expression of the dephasing rate, computed in Ref. [19]. 



3 Anderson model 

Let us next discuss the Anderson model, defined by the Hamiltonian, 

H = J2 e( P K CT a pCT + e d J2 d\d a + Ud\d } d\d x + V £ (c+A + h.c.) . 

per cr CT,p 

Here d a denotes a local ci-level's annihilation operator, U is the on-site Coulomb 
repulsion, and the conduction band and the local electronic level are hybridized 
by V. The Anderson model is the most elementary Hamiltonian that describes 
local moment formation on a localized ci-orbital, and in fact, the Kondo Hamil- 
tonian can be obtained from it in the limit A U + e^, \ed\, with A = irgV 2 
the width of the resonance [29]. 

As first discussed by Langreth [27], the T-matrix for the Anderson model can 
be related to the (i-level's Green's function as [9,29] 

Im{T»} = W 2 Qd ^{uj) , Re{T a (u)} = -r V 2 f dJ . (16) 

J uj — UJ 

Here r = sgn u>, and Qd,a(w) is the spectral function of the <i-Fermion, which 
we have computed using NRG [21]. 
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Fig. 4. Elastic, inelastic and total scattering cross section for the asymmetric An- 
derson model. The low-energy part of the curves is very similar to the one obtained 
for the Kondo model. 

The full frequency-dependence of the various scattering cross-sections obtained 
for the asymmetrical Anderson model with interaction strength U / A = 5.1 is 
shown in Fig. 4. At this value of U one can already observe the Hubbard side- 
peaks in the total scattering cross section at energies to ~ and uj « U + e^, 
and a distinct Kondo resonance appears at uj w too. The scattering rates 
in the region uj w are strikingly similar to the ones we obtained for the 
Kondo model, and remarkably, both the quasi-linear regime of (Xi ne i and the 
plateau are already present for these moderate values of U/A. This is not 
very surprising since, as stated before, the Kondo model is just the effective 
model of the Anderson model in the limit of large U/A and uj U. For 
even larger values of U/A and intermediate energies, T K <^ uj <^ U, the 
elastic and inelastic contributions follow very nicely the asymptotic behavior 
found for the Kondo model, and scale as ~ 1/\u\uj/T k ) and ~ 1/ ln 2 (uj/T K ), 
respectively. New features compared to the Kondo model are the Hubbard 
peaks that correspond almost entirely to inelastic scattering. 



4 Inelastic scattering in the two-channel Kondo model 

So far, we considered scattering from Fermi liquid models only. Let us now 
discuss the two-channel Kondo model, the prototype of all non-Fermi liquid 
impurity models [28]. This is defined by a Hamiltonian similar to Eq. (10) 
excepting that now there is two 'channels' of conduction electrons, a = 1, 2 
that are coupled to the impurity spin with couplings J a , 

H = J2^P fl UV,a+ -^^Y, a i,a°™' a p',a • (17) 

0=1,2 p,cr 0=1,2 L p.p' 

an* 
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Fig. 5. Left: Renormalization group flow of the eigenvalue of the single particle 
S'-matrix for the two-channel Kondo model. Right: Energy-dependence of elastic and 
inelastic scattering rates for the two-channel Kondo model in units of <7o = 4-ir/p'p, 
at T = 0. 

In the channel-symmetric case, J\ = J2 the two conduction electron channels 
compete to screen the impurity spin independently, which is therefore never 
completely screened. This competition leads to the formation of a strongly 
correlated state which cannot be described by Nozieres' Fermi liquid theory, 
and is characterized by a non-zero residual entropy, the logarithmic diver- 
gence of the impurity susceptibility, and the power law behavior of transport 
properties with fractional exponents [28] . Any infinitesimal asymmetry in the 
couplings A = {J 1 — J2)/{J\ + J 2) leads to the appearance of another low- 
temperature energy scale T* oc A 2 Tk/ at which the system crosses over to a 
Fermi liquid behavior: Electrons being more strongly coupled to the impurity 
form a usual Kondo singlet with the impurity spin, while the other electron 
channel becomes completely decoupled from the spin. 

For Ji = J2, no Fermi-liquid relations are available. There exists, however, 
an exact theorem due to Maldacena and Ludwig, according to which, at the 
two-channel Kondo fixed point, the single-particle elements of the S'-matrix 
identically vanish for uj — > 0: s 2 ck(u — ► 0) = [32]. As a consequence, 
t-icK^y 1 — 0) — —i ■ This relation leads to the surprising result that exactly 
half of the scattering is inelastic at the Fermi energy, while the other half of 
it is inelastic: 

°H<(" = 0) = ofV = 0) = ogV = 0)/2 . (18) 



This counter-intuitive result can be understood as follows: The vanishing of the 
single particle S'-matrix indicates that an incoming electron cannot be detected 
as one electron after the scattering event, and it "decays" into infinitely many 
electron-hole pairs. To get such a "decay", however, the scattering process 
must have an elastic part which interferes destructively with the unscattered 
direct wave, and cancels exactly the outgoing single particle amplitude in the 
s-channel. 
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The evolution of s(u) and the inelastic scattering rates for the two-channel 
Kondo model are shown in Fig. 5 as a function of the energy of the incoming 
particle. In the channel-symmetric case inelastic processes are allowed even at 
u — 0, which is a clear signature of the non- Fermi liquid behavior. The non- 
Fermi liquid nature is also reflected in the ~ -Juj singularity of the scattering 
cross sections at uo = 0. 

For A > the total scattering rate approaches the unitary limit in channel 
"1" below the Fermi liquid scale T*, while it goes to for A < 0. In both 
cases, the inelastic scattering freezes out, o" inc i(u;) shows a dip below T*, and 
it ultimately scales to as 

<7inel(u;) OC UJ 2 /T* 2 . 

Remarkably, the inelastic scattering cross-section is very similar for A > and 
A < 0, while the scattering contributions are dramatically different in these 
two cases (see Ref. [21]). 



5 Conclusions 

In the present paper we reviewed a recent theory of inelastic scattering from 
quantum impurities, and applied it to the single and two-channel Kondo mod- 
els, and the Anderson model. 

We showed that in the Kondo model and in the local moment regime of the 
Anderson model a broad plateau appears in the T = temperature energy- 
dependent inelastic scattering rate above Tk, while a quasi-linear regime 
emerges below Tk, in the T = temperature, both in excellent agreement 
with recent experimental observations. We also computed the universal flow 
of the eigenvalues of the ^-matrix, which is a useful quantity to classify various 
types of impurity states [31]. 

As an example of a non- Fermi liquid, we also discussed scattering from the two- 
channel Kondo model, where half of the scattering remains inelastic even at the 
Fermi energy, and correspondingly, the eigenvalue of the S'-matrix vanishes. 
This fragile non- Fermi liquid state is, however, destroyed, once a small channel- 
symmetry breaking is applied, and then the scattering becomes elastic below 
a Fermi-liquid scale, T* [21]. 

Our results have been extended to T ^ temperatures by Micklitz et al. 
in Ref. [19], who showed that the dephasing time can indeed be related to 
the inelastic cross-section computed here for T = temperature. The theory 
presented here has also been applied by Koller et al. to the S = 1 Kondo 
model [35]. This under-screened model represents a singular Fermi liquid [31], 
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where the scattering is elastic at the Fermi energy, however, it scales to only 
logarithmically, in contrast to a Fermi liquid. 

Let us close our conclusions with an important remark. In a real experiment, 
the external electromagnetic field couples with a minimal coupling to the con- 
duction electrons. As a consequence, the Kubo formula is formulated in terms 
of the conduction electron current operators, and the relevant quantity to de- 
termine dephasing is thus the inelastic scattering rate of electrons. This is 
what we compute here and that has been computed in Ref. [19]. Quasiparti- 
cles are, on the other hand, typically not minimally coupled to the guage field, 
since they are usually complicated objects in terms of conduction electrons. 
If one defines quasiparticles as stable elementary excitations of the vacuum, 
as Nozieres did [23], or as they appear in Bethe Ansatz, then, by definition, 
these quasiparticles do not decay at all at T = and scatter only elastically 
[23]. However, excepting for u = 0, a real conduction electron is composed 
of many such stable quasiparticles, and it already decays inelastically even at 
T = temperature. In the Kondo model, at the Fermi energy quasiparticle 
states are just phase shifted conduction electron states, however, the connec- 
tion between quasiparticles and conduction electrons is not trivial for any 
finite energy. Therefore, if one considers inelastic scattering at a finite energy, 
one must precisely specify how finite energy quasiparticle states are defined, 
how they couple to a guage field, and how a finite energy electronic state is 
decomposed in terms of these quasiparticles,. In the present framework, we 
avoid this difficulty by formulating the problem in terms of electrons. 

We are indebted to L. Saminadayar, C. Bauerle, J.J. Lin, and A. Rosch for 
valuable discussions. This research has been supported by Hungarian grants 
Nos. NF061726, D048665, T046303 and T048782. L.B. acknowledges the fi- 
nancial support of the Bolyai Foundation. 
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